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ABSTRACT 

This paper examines the generation of seed magnetic fields due to the growth of cosmological 
perturbations. In the radiation era, different rates of scattering from photons induce local differences 
in the ion and electron density and velocity fields. The currents due to the relative motion of these 
fluids generate magnetic fields on all cosmological scales, peaking at a magnitude of 0(10^^"' G) at 
the epoch of recombination. Magnetic fields generated in this manner provide a promising candidate 
for the seeds of magnetic fields presently observed on galactic and extra-galactic scales. 
Subject headings: cosmology : theory — early universe — large-scale structure of universe — galaxies: 
magnetic fields 



1. INTRODUCTION 

The presence of magnetic fields on galactic and ex- 
tragalactic scales is a major unsolved puzzle in mod- 
ern astrophysics. Although the observational evidence 
supporting the existence of magnetic fields in large-scale 
structure is overwhelming, there is no consensus as to 
their origins. The standard paradigm for the creation of 
these fields is the dynamo mechanism, in which a small 
initial seed field is amplified by turbulence (the a effect) 
and/or by differential rotation (the oj effect). In princi- 
ple, once a seed field is in place, it should be possible to 
follow its evolution and amplification from the collapse 
of structure and the operation of any relevant dynamos. 

In recent papers, it has been argued that magnetic 
fields are generated in the early universe due to induced 
vorticity (Matarrese et al. 2005; Gopal & Sethi 2005; 
Takahashi et al. 2005). In this paper, we put forward a 
new mechanism for the generation of seed fields from cos- 
mological perturbations. We show that the evolution of 
cosmological perturbations in the pre-recombination era 
produces charge separations and currents on all scales, 
both of which contribute to magnetic fields, in addition 
to and independent of any fields generated via vortic- 
ity. These seed fields persist until the onset of significant 
gravitational collapse, at which point field amplification 
and dynamo processes may magnify such seeds to the 
0(^0) fields observed today. 

The generation of magnetic fields from charge separa- 
tions and currents in the early universe is a necessary 
consequence of structure formation. This paper calcu- 
lates the magnitude of these seed fields, and it is shown 
that these seed fields may be sufficiently strong to ac- 
count for all of the observed magnetic fields in large-scale 
structures. The layout of the paper is as follows: Sec- 
tion 2 gives an overview of the observational evidence for 
magnetic fields along with a brief theoretical picture of 
their generation. Section 2 also contains an explanation 
of the novel idea that the early stages of structure for- 
mation in a perturbed universe generate magnetic fields. 
Section 3 presents a detailed treatment of cosmological 
perturbations, with a specific view towards the creation 
and evolution of local charge separations and currents. 
Section 4 explicitly calculates the seed magnetic fields 



that arise via this mechanism as a function of scale and 
epoch. Finally, section 5 compares the results of this 
mechanism with competing theories. It also discusses 
avenues for future investigation of this topic, including 
possible observational signatures of the fields that would 
arise from this mechanism. 

2. MAGNETIC FIELDS: BACKGROUND 

Magnetic fields with strength ^ jjlG are seen in all 
gravitationally bound or collapsing structures in which 
the appropriate observations have been made (Widrow 
2002). The four major methods used to study astrophys- 
ical magnetic fields are synchrotron radiation, Faraday 
rotation, Zeeman splitting, and polarization of starlight. 
These observational techniques are detailed in depth in 
Ruzmaikin et al. (1988), with Faraday rotation often 
proving the most fruitful of the above methods. 

Magnetic fields have been found in many different 
types of galaxies, in rich clusters, and in galaxies at high 
redshifts. Spiral galaxies, including our own, appear to 
have relatively large magnetic fields of 0{\Q jiG) on the 
scale of the galaxy (Fitt & Alexander 1993), with some 
(such as M82) containing fields up to ~ 50 ^G (Klein et 
al. 1988). Elliptical and irregular galaxies possess strong 
evidence for magnetic fields (of order ~ /j,G) as well 
(Moss & Shukurov 1996), although they are much more 
difficult to observe due to the paucity of free electrons in 
these classes of galaxies. Coherence scales for magnetic 
fields in these galaxies, as opposed to spirals, are much 
smaller than the scale of the galaxy, of order 100 to f 000 
pc. Furthermore, galaxies at moderate (z ~ 0.4) and 
high redshifts [z > 2) have been observed to require sig- 
nificant {^ fxG) magnetic fields to explain their observed 
Faraday rotations (Kronberg et al. 1992; Athreya et al. 
1998). Magnetic fields are also observed in structures 
larger than individual galaxies. The three main types of 
galaxy clusters are those with cooling flows, those with 
radio-halos, and those devoid of both. Galaxy clusters 
with cooling flows are observed to have fields of 0.2 to 
3/iG (Taylor et al. 1994), the Coma cluster (a prime ex- 
ample of a radio-halo cluster) is observed to have a field 
strength r^ 2.5 ^G (Kim et al. 1990), while clusters se- 
lected to have neither cooling flows nor radio halos still 



exhibit indications of strong {0.1-1 iiG) fields (Clarke et 
al. 2001). There even exists evidence for magnetic fields 
on extracluster scales. An excess of Faraday rotation is 
observed for galaxies lying along the filament between 
the Coma cluster and the cluster Abcll 1367, consistent 
with an intercluster magnetic field of 0.2-0.6 /iC (Kim 
et al. 1989). On the largest cosmological scales, there 
exist only upper limits on magnetic fields, arising from 
observations of the cosmic microwave background (Bar- 
row et al. 1997; Yamazaki et al. 2005, 2006) and from 
nucleosynthesis (Cheng et al. 1994), setting limits that 
on scales > lOMpc, field strengths are < 10~^ G. 

Observational evidence for magnetic fields is found in 
galaxies of all types and in galaxy clusters, both locally 
and at high redshifts, wherever the appropriate observa- 
tions can be made. A review of observational results can 
be found in Vallee (1997). The theoretical picture of the 
creation of these fields, however, is incomplete. Fields of 
strength ~ /iG can be explained by the magnification of 
an initial, small seed field on galactic (or larger) scales 
by the dynamo mechanism (Parker 1971; Vainshtein & 
Ruzmaikin 1971, 1972). A protogalaxy (or protocluster) 
containing a magnetic field can have its field strength 
increased by many orders of magnitude through gravita- 
tional collapse (Lesch & Chiba 1995; Howard & Kulsrud 
1997), and can also be further amplified via various dy- 
namos. Dynamos which can amplify a small seed field 
into the large fields observed today involve helical turbu- 
lence (a) and/or differential rotation {lo). Various types 
of these dynamos include the mean-field dynamo (Steen- 
beck et al. 1966; MoflPatt 1978; Ruzmaikin et al. 1988), 
the fluctuation dynamo (Kazantsev et al. 1985; Moss 
& Shukurov 1996), and merger-driven dynamos (Tribble 
1993), among others. However, the dynamo mechanism 
does not explain the origin of such seed fields. 

While the initial seeds that grow into magnetic fields 
are anticipated to be small, they must still come from 
somewhere (Zel'dovich & Novikov 1983), as their exis- 
tence is not explained by the dynamo mechanism alone. 
There are many mechanisms that can produce small- 
strength magnetic fields on astrophysically interesting 
scales, either through astrophysical or exotic processes 
(see Widrow 2002, for a detailed review). Exotic pro- 
cesses generally rely on new physics in the early universe, 
such as a first-order QCD phase transition (Hogan 1983; 
Quashnock et al. 1989), a first order electroweak phase 
transition (Baym et al. 1996; Sigl et al. 1997), broken 
conformal invariance during inflation (Turner & Widrow 
1988; Ashoorioon & Mann 2005), speciflc inflaton po- 
tentials (Ratra 1992), or the presence of charged scalars 
during inflation (Calzetta et al. 1998; Kandus et al. 2000; 
Davis et al. 2001). Astrophysical mechanisms, in con- 
trast, are generally better grounded in known physics, al- 
though they have difficulty generating sufficiently strong 
fields on sufficiently large scales. The difference in mo- 
bility between electrons and ions admits, under appro- 
priate circumstances, the creation of large-scale currents 
and magnetic fields. Seed magnetic fields arising from 
this difference in mobility between electrons and ions 
may be astrophysically generated via a variety of phe- 
nomena, including radiation-era vorticity (Harrison 1970, 
1973), vorticity due to gas-dynamics in ionized plasma 
(Biermann 1950; Pudritz & Silk 1989; Subramanian et 
al. 1994; Kulsrud et al. 1997; Gncdin et al. 2000), stars 



(Syrovatskii 1970), or from active galactic nuclei (Hoyle 
1969). Although there are many candidates for pro- 
ducing the seed magnetic fields required by the dynamo 
mechanism, none has emerged as a definitive solution to 
the puzzle of explaining their origins. 

The mechanism proposed in this paper is that seed 
magnetic fields are generated by the scattering of photons 
with charged particles during the radiation era. Unlike 
the mechanism of Harrison (1970, 1973), which is disfa- 
vored (Rees 1987) due to its requirement of substantial 
primordial vorticity (although see Matarrese et al. 2005; 
Takahashi et al. 2005, for an argument that some vortic- 
ity is necessary), the fields of interest here are generated 
by the earliest stages of structure formation, requiring no 
new physics. Ions (taken to be protons, for simplicity) 
and electrons are treated as separate fiuids, with oppo- 
site charges but significantly different masses. The mass- 
weighted sums of their density and velocity fields will 
determine the evolution of the net baryonic component 
of the universe, and should agree with previous treat- 
ments (e.g. Ma & Bertschinger 1995). The differences 
of their density and velocity fields, however, provides a 
description of a local charge separation and a local cur- 
rent density, both of which contribute to magnetic fields. 
Since cosmological perturbations, which serve as seeds 
for structure formation, exist on all scales, it is expected 
that seed magnetic fields will be generated on all scales 
by this mechanism. As the magnetic fields generated via 
this mechanism are in place at very early times, they 
possess the advantage over other astrophysical mecha- 
nisms that these fields will have optimally large amounts 
of time to be affected by amplification effects and var- 
ious dynamos. The remainder of this paper focuses on 
calculating the magnitude of the magnetic fields gener- 
ated by this process and discussing their cosmological 
ramifications. 

3. COSMOLOGICAL PERTURBATIONS 

Although the early universe is isotropic and homoge- 
neous to a few parts in lO'"' (Bunn & White 1997), it is 
these small density inhomogeneities, predicted by infia- 
tion to occur on all scales (Guth & Pi 1982), that lead 
to all of the structure observed in the universe today. 
As it is the early epoch of structure formation that is of 
interest for the creation of magnetic fields, we calculate 
the evolution of inhomogeneities in the linear regime of 
structure formation. The most sophisticated treatment 
of cosmological perturbations in the linear regime to date 
is that of Ma & Bertschinger (1995), which provides evo- 
lution equations for an inhomogeneous universe contain- 
ing a cosmological constant, dark matter, baryons, pho- 
tons, and neutrinos. This section extends their treatment 
to encompass separate proton and electron components. 
The mass- weighted sum of protons and electrons will re- 
cover the usual baryon component, whereas the differ- 
ence of the density fields represents a charge separation, 
and the difference of the velocities is a net current. 

The dynamics of any cosmological fluid can be ob- 
tained, in general, from the linear Einstein equations 
(see Peebles & Yu 1970; Silk & Wilson 1980; Wilson & 
Silk 1981, for earlier treatments). Although the choice of 
gauge does not impact the results, the Conformal New- 
tonian gauge leads to the most straightforward calcula- 
tions. When tensor modes are unimportant, the metric 



47r 



is given by 

ds^ ^ a^(T)[-(l + 2ij)dT^ + (1 - 2(l))dx'dx^], (1) 

where ip :^ (j) when gravitational fields are weak. The 
linearized Einstein equations are then 

fc20 + 3 -(</.+ - ^) =. AirGa^STl, 
a a 

fc2 (</) + - V) - 47rGa2(p + P)e, 
a 

4> + -{i^ + 24>) + i—- ^)^ + ^{^ - ^) 

a a a^ 6 

k'^{(t, - V') = 127rGa^(p + P)ct, 

where a is the shear stress, which is negligible for non- 
relativistic matter but important for photons and neutri- 
nos. A cosmological fluid that is either uncoupled to the 
other fluids or mass-averaged among uncoupled fluids in 
the early universe obeys 

5 = -(l + w){e- 30) - 3-(c2 - w)5, 
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where 5 is defined as the local density relative to its spa- 
tial average {S = Sp/p); 9 = ik^Vj, where Vj is the local 
peculiar velocity; and Cg is the sound speed of the fluid. 
For individual components with additional inter- 
component interactions, equation (3) must be augmented 
to include these interactions. Examples of such interac- 
tions include momentum transfer between photons and 
charged particles and Coulomb interactions between pro- 
tons and electrons. For protons, electrons, and cold dark 
matter (CDM), an equation of state w = is assumed, 
and for radiation and neutrinos, w ~ ■^. The master 
equations for each component of interest are computed 
explicitly in subsections 3.1-3.5. The evolution equa- 
tions which appear in this section for baryons, photons, 
light neutrinos, and cold dark matter are derived in much 
greater detail, including photon polarization, hydrogen 
and helium recombination, in both synchronous and con- 
formal Newtonian gauges, in Ma & Bertschinger (1995). 

3.1. Cold Dark Matter 

The cold dark matter component (denoted by the sub- 
script c) is collisionless and pressureless; its evolution can 
be simply read off from equation (3) with w = c^ = a = 
to be 
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k^ip. 



(4) 



Any cold (i.e. nonrelativistic), collisionless component 
will behave according to the dynamics given by equation 

(4). 

3.2. Light Neutrinos 

For massless (or nearly massless) particles, pressure 
is non-negligible, and the appropriate equation of state 
has w ~ \^. Additionally, the shear term {a) may be 
important as well. The only accurate way to compute 



the evolution of such a component of the universe is by 
integration of the Boltzmann Equation, which is given 
for light neutrinos (denoted by subscript v) by 



dr 



+ ik{k ■ fi)J-u = A[(j) — ik{k ■ n)?/;]. 



(5) 



in Fourier space. 

The approximation that neutrinos are massless and un- 
coupled is very good from an age of the universe of ap- 
proximately i ~ 1 s until the epoch of recombination. 
The evolution equations for light neutrinos are then 



<5,. = --e,. + #, 



e^^k^{-5^-cT„+^), 



T, 



w = ^^-pY[^^i.(i-i) - (? + l).F^((+i)], (6) 

where Oi, is related to J-y by 2a i, = Jv2, and the index 
/ governs the final equation for I > 2. Tui is defined 
by the expansion of the perturbations in the distribution 
function, J-^, 

oo 

J-, = ^(-z)'(2Z + l)J-,,(fc,T)P,(fc-n), (7) 

;=o 

where Pi{k ■ fi) are the Lcgcndre polynomials. Equations 
(5-7) are valid for any non-coUisional species behaving 
as radiation. 

3.3. Photons 

Photons (denoted by subscript 7), although similar 
to light neutrinos, evolve differently due to their large 
coupling to charged particles. The differential scattering 
cross-section of photons with electrons, is given by 



da 3(Tt 



cos^e). 



(8) 



where ctt is the Thomson cross-section (Jackson 1998). 
Photons also scatter with protons, but with a cross- 
section suppressed by a factor of the squared ratio of 
electron to proton mass, ml/rai « 3.0 x 10~^. 

The polarization-summed phase-space distribution for 
photons, ^-y, is the same as the distribution function for 
neutrinos (see equation 7) , with a non-zero difference be- 
tween the two linear polarization components denoted by 
Q^. The linearized collision operators for Thomson scat- 
tering (Bond & Efstathiou 1984, 1987; Kosowsky 1996; 
Ma & Bertschinger 1995) yield the set of master equa- 
tions for photons. 
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where JF^o — 5^, -^71 — 4:9^/3k, T^i 
dices I and to are valid for Z > 3 and rri > 0, and 
the subscript h denotes the baryonic component, which 
is the mass-weighted sum of the electrons and protons. 
Electron-photon scattering is so dominant over proton- 
photon scattering as to render the latter negligible, but 
the electron-proton coupling (via electromagnetism) is 
sufficiently strong that, to leading order, those two fluids 
move in kinetic equilibrium. 

3.4. Baryons 

The net behavior of the baryonic component can be 
derived from combining the mass-weighted contributions 
of the proton and electron fluids. Both protons and 
electrons contain all of the terms present in the CDM 
equations (cf. Section 3.1), but additionally contain im- 
portant terms arising from the Coulomb interaction and 
from Thomson scattering. The coupling of the Coulomb 
interaction to density inhomogeneities can be calculated 
through a combination of the electromagnetic Poisson 
equation, 

V2$ = -V • i? = -47rpc, (10) 

where pc is the electric charge density, and the Euler 
equation, 

i^ + i(.^.v).--iv0-^iva> + c, (11) 

a at a a ma 

with q/m as the charge-to-mass ratio of the particle in 
question and C the collision operator. The Coulomb con- 
tribution appears as Airea? (np—ne)cii / m,i in the evolution 

equation for 9i , where i denotes a species of particle with 
a mass mi and charge q^. The evolution equations are 
therefore 



5e = -t 



30, 



Oe = ee+cik^Se+k^^ 

a 

+ Te{e^-ee) (rip-ne), (12) 

TOe 

for electrons (denoted by subscript e), and 

5p = -0p + 3(j), 



9p^ Op + cj.k^5p + k^ip 



r. 



Attc^o^ 
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-{np-He), (13) 



for protons, denoted by subscript p, where F character- 
izes the rate of momentum transfer due to photon scat- 
tering with charged particles. The damping coefhcient 
for electrons (Fg) is given by 



Fe = 



_Ap^ne<7Ta 



3/Oe 



(14) 



and the analogous quantity for protons is smaller by a 
factor Fp/Fe = [me/mpf 2± 1.6 x lO-^". Note the dif- 
ference in the sign of the final terms in the equations for 
^e and 9p, which will prove important below. 



From equations (12) and (13) for electrons and protons 
the dominant gravitational and electromagnetic combi- 
nations can be constructed separately. The remainder 
of this subsection details the evolution of baryons in the 
linear regime of a perturbed universe. Baryonic matter 
can be treated as the combination of electrons and pro- 
tons, thus the mass weighted sum of proton and electron 
overdensities gives rise to the baryonic perturbations. 
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where mi, 
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TTig. The evolution of baryonic matter 



follows from a mass-weighted combination of the equa- 
tions for electrons (12) and protons (13). To the ex- 
tent that electrons and protons move together (the tight- 
coupling approximation) and the difference Up—Ue is neg- 
ligibly small (for baryonic evolution this is an excellent 
approximation; see Section 3.5), the evolution equations 
for baryonic matter are 



-Ob + 3(1), 

--9b + clkHb + k^tl^ + Tb{9-, 



(16) 



where F;, = (TO-pFp -I- meTe)/mb ~ Tem,e/mb. In the 
tight coupling approximation, the baryon-photon cou- 
pling term in equation (16) is driven by the electron- 
photon interaction, as has been shown by Harrison (1970) 
and subsequent authors. The equations in (16) are iden- 
tical to the equations for the evolution of baryonic inho- 
mogeneities derived in Ma & Bertschinger (1995). 

3.5. Charge Separation 

From equations (12) and (13), a charge difference com- 
ponent as well as a sum component can be obtained. As 
the limits on a net electric charge asymmetry in the uni- 
verse are very strict (Orito & Yoshimura 1985; Masso 
& Rota 2002; Caprini & Ferreira 2005), we expect that 
any differences in densities and/or velocities of protons 
and electrons will not be strong enough to significantly 
alter the evolution of the other components of the uni- 
verse. It is therefore expected that the analysis of bary- 
onic matter, detailed in Section 3.4, will be unaffected 
by differences in proton and electron densities and veloc- 
ities. However, even a small separation may still have 
important cosmological implications. 

The charge difference component (denoted by sub- 
script q) is the difference between the proton and elec- 
tron components, such that Sq = Sp — dg and 9q = 9p — 9e. 
The gravitational potential does not affect the evolution 
of these quantities; gravity acts equivalently on electrons 
and protons. However, if Sq and 9q are nonzero, result- 
ing electromagnetic fields will act differently on oppo- 
sitely charged particles. The master equations for the 
charge- asymmetric component, obtained from equations 
(12) and (13) are 



Sq = -9q 



o2„2 



9q^ Oq + clPSq + 6q 

a TO-e 



Fe(07 — 9b+ 9q), 



(17) 



where the approximations Fp <C Fe and to,;, ~ rrip have 
been utilized where applicable. The term Aimee^a^Sq/me 



in equation (17) arises from the Coulomb force acting on 
charged particles, while the final term, re{9^ — Ob + Oq), 
arises from the difference in Thomson scattering between 
protons and electrons. This final term provides a source 
for the generation of a charge separation independent of 
and in addition to any initial charge asymmetry, and will 
create a local charge asymmetry even when there is none 
initially. In the evaluation of equation (17), the electro- 
magnetic terms dominate the cosmological terms, such 
that an excellent approximation in the pre-recombination 
universe is 



4TTnf,e^a^ 



(18) 



For some purposes, it is useful to express the set of 
equations found in equation (17) as a single ordinary dif- 
ferential equation. This can be accomplished by setting 
9q — —Sq, and again by neglecting the unimportant cos- 
mological terms Oqh/a and c^k^Sq. Many of the coeffi- 
cients in equation (17) are functions of a, but the deriva- 
tives in equation (17) are with respect to conformal time, 
T. A change of variables can be performed, using the re- 
lation that 
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Abh^c 



327r3G(fcr)4 

N 



-Na'^to 



1/2 



72.2, 



(19) 



in the radiation era, where to is the age of the universe 
today and t = ar, to express all derivatives as derivatives 
with respect to a, denoted by primes (instead of dots). 
The evolution of 5q can be tracked by evolving equation 
(20) below. 
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(20) 



where the subscript denotes the present value of a quan- 
tity. This is simply the equation of a damped harmonic 
oscillator, with coefficients that change slowly with time 
compared to damping or oscillation times. The behav- 
ior can be characterized as overdamped at the earliest 
times, critically damped when a « 2.3 x 10^^, and un- 
derdamped at late times. 

Of all the terms in equation (20), only 9^, 9b, and Sq 
(and derivatives) are functions of a; all other quantities 
are constant coefficients. Although there does not exist a 
simple analytic form for [9^~9b) in general, at sufficiently 
early times there exists the simple approximation 

9^-9h^ 6.0 X 10^^ fc" 



'^^a\ 



(21) 



valid when the following condition is met 
a < 



10"^ for k < 0.1Mpc"\ 



< 10" 



IMpc"^ 



for k > 0.1 Mpc" 



Equation (21) is an approximation for a flat 
ACDM cosmology with cosmological parameters 
Hq = 71kms-iMpc"\ f]„ = 0.27, Qb = 0.044, and a 
Helium-4 mass fraction oi Y = 0.248. These parameters 
are used in all subsequent analyses for the calculation of 
cosmological quantities. 



The approximation in equation (21) breaks down at 
later times. When this occurs, numerical methods must 
be used to obtain the quantity {9^ — Ob). The software 
package CDSMICS (Bertschinger 1995) is ideal for per- 
forming this computation, as it performs the numerical 
evolution of equations (4,6,9, and 16) concurrently. Com- 
putational results for the quantities 9j and 9b are given 
by COSMICS, which are valid at all times in the linear 
regime of structure formation. It is found that when the 
approximation in equation (21) breaks down, the quan- 
tity {9-y — 9b) grows more slowly initially, and proceeds to 
oscillate at a roughly constant amplitude at later times. 
These oscillations in the quantity {9^ — Ob) are closely 
related to the acoustic oscillations between baryons and 
photons observed in the cosmic microwave background 
(Bennett et al. 2003). 

Numerical integration of equation (20) can be accom- 
plished in various ways (see Press et al. 1992, for exam- 
ples). At sufficiently late times (when a ^ 2.3 x 10~^), 
numerical results indicate that the quasi-equilibrium so- 
lution 

tlA fl^\(^g 0^^^ (22) 

obtained by neglecting the first two terms in equation 
(20), is an excellent approximation. With 9j and 9b given 
by COSMICS in units of Mpc""'^, the prefactor in equation 
(22) can be written as 
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(23) 



where re = e^ /meC^ is the classical electron radius. The 
quantity 9q then follows directly from equation (17) to 
be 

crrmb f Py,o^ ( }_ 

i2 



3 7r e^ 



P6,0 



'1 



0- 



(24) 



The solutions in equations (22) and (24) remain 
valid up to the epoch of recombination [z ~ 1089). 




Fig. 1. — Fourier mode amplitudes of the charge asymmetry 
Sq plotted as a function of redshift z for comoving scales of (from 
top to bottom) k = 10, 1, 0.1, and 0.01 Mpc~^. The amplitude 
5q rises as ~ a^ initially, then ceases to grow when the scale of 
interest enters the horizon, and oscillates at an amplitude which 
first continues to rise slowly, then falls, eventually matching on 
to the equilibrium solution &q oc 6;,. Solid lines indicate positive 



The results of numerically integrating the equations 
for 5q and 9q on various length scales through recom- 
bination using the code COSMICS (Bertschinger 1995) 
are presented as a function of rcdshift z in Figure 1 for 
k = 10, 1, 0.1, and 0.01 Mpc"^ The amplitude 5q rises 
as ^ a^ initially, then ceases to grow when the scale of 
interest enters the horizon, and oscillates at an amplitude 
which first continues to rise slowly, then falls, eventually 
matching on to the equilibrium solution Sq oc 9b- 

The magnitude of the charge separation is small (of or- 
der 10"^"^ in Figure 1) because the constant in equation 
(23) is small compared to the Hubble length at decou- 
pling. In the conformal Newtonian gauge the mode am- 
plitudes are defined to be of order 1 outside the horizon, 
and need to be multiplied by the COBE normalization 
Sh = 1.95 X lO"'^ of Bunn & White (1997) to obtain the 
physical quantities of magnetic field strength and mag- 
netic spectral density, as in Section 4. 

Figure 2 shows the spectral density of the charge 
asymmetry, iirk^ Pq{k) / {2tt)^ , plotted as a function of 
k. For comparison, the spectral density of baryonic mat- 
ter, 47rfc^P6(fc)/(27r)'^, is also shown. Note the identical 
scales, fe, at which Pq{k) and Pb{k) change signs. 




Fig. 2. — Spectral density of the charge asymmetry, 
47rk^Pq{k)/{2TT)^ , at the epoch of decoupling, z = 1089, plotted 
as a function of wavenumber k (lower graph) . Also shown, at the 
same epoch, is the spectral density of baryons, 47rfc'^Pi,(fc)/(27r)^, 
(upper graph). 

The results thus far are accurate up through the epoch 
of recombination. At this point, however, the universe 
transitions from a state with ionization fraction Xe ~ 1 to 
a state where the ionization fraction is small, Xe ~ 10"'^ 
(Peebles 1968), and the universe is transparent to radia- 
tion. The calculations of Matarrese et al. (2005) suggest 
that Silk damping affects only scales k > lMpc~^. In 
the absence of interactions with photons, a charge sepa- 
ration would continue to evolve as 

dq = dq = --IhH K a dq, (25) 

where K^ is the ratio of the electric to gravitational 
forces, 

K^ = = 2.26 X 10^^. 

Grripme 



A residual charge separation without additional interac- 
tions would oscillate (plasma oscillations) with proper 
frequency to « KH, with a further slow adiabatic decay 
of the oscillation amplitude. However, there are many 
other effects that begin to become important after re- 
combination, including gravitational collapse, dynamo 
effects, and continued electron-photon scattering, and we 
do not attempt to compute the subsequent evolution in 
detail in this paper. 

4. MAGNETIC FIELDS 

With the results derived in Section 3.5 for 6q and 6q, 
values for the local current densities and local charge 
separations can be obtained at any time in the pre- 
recombination universe on all scales. Both 6q and 9q will 
contribute to magnetic fields, as currents create mag- 
netic fields directly, and the bulk motion of a region of 
net charge will also produce a magnetic field. For each 
comoving distance scale (given by the value of k) and 
each epoch (determined by the scale factor a) of the uni- 
verse, there will be a magnetic field amplitude associated 
with that scale. This field may serve as the seed for the 
large-scale magnetic fields observed today. 

An expression for magnetic fields can be derived from 
the currents arising from the relative motion of the pro- 
tons and electrons in the universe. Magnetic fields can 
be derived from Maxwell's equations 



V-B = 0, V X B = An J 



% '-) 



with the current density J given by 

n^eVe^nee[5qVb + {1 + 5b)Vq\, (27) 



fj '''O ^ ^T) 



where Vq = Vp — v,,, and the displacement current is ne- 
glected. 

In the Fourier domain, the curl of Maxwell's equation 
for V X i? gives a direct expression for the magnetic field 
as a convolution 



m=- 



AnUgQe f (Pk' k x k' 



a2|^|2 J (2^)3 |fc/|2 

[9bik') Sq{k - k') + 9q{k') Sbik - k% (28) 

from which the second moment of the magnetic field B{k) 
is 



{B,{ki)Bjik2)) = i27rfSuiki + fcs) ^^ ^^(fc), 



(29) 



where Su is the Dirac delta function and Psik) is the 
magnetic field power spectrum. Note that the direc- 
tion parallel to k does not contribute to magnetic fields. 
Therefore, the direction perpendicular to k is projected 
out in equation (29). The power spectrum, Psik), is 
then given by 



Psik)-^ 




2 A' 



(27r)3 ' 
,oMk'\)P5.sMk-k'\) 



PeM\k'\)PeM\k~k'\) 



\k-k'\^ 
^ PeM\k'\)Ps,5Mk-k'\) 



l/c'P 



1 



|A:-fc'|2 
2 



^M,(l^'l)^M,(|fc-^'l) 



PeM\k'\)PsM\k-k'\) 



Fig. 3. — Spectral density of the magnetic field 47rfe''Ps(A;)/(27r)^ 
(in G^) generated by cosmological perturbations on a given comov- 
ing scale k (in Mpc"'^) at the epoch of recombination z ~ 1089. 
The dashed lines show the three separate contributions in equation 
(30) (where the term containing Pg g is short-dashed; Pe^^e), is 
long-dashed; and Pg q^ is short/long-dashed), and the solid (top- 
most) line shows the total magnetic spectral density. The peak 
value corresponds to a magnetic field strength B ~ 10"-^^— 10~'^* G. 
Note also that the oscillations due to Silk damping become unim- 
portant for k > 0.5 Mpc^^. 



|/c-fc'|2 



Ai\k'\)PB.sA\k~k'\) , (30) 



where A' is the angle between the vectors k and k'. The 
second moment of two general quantities, (j) and ip, is 
expressed as the power Ptp^{k), defined by 



(^(fci)V^(fc2)) = (2^)^<5D(fci + fca) P^i^ik), (31) 



where k = \ki\ = \k2\. All terms in equation (30) make 
comparable contributions to the total. Note that al- 
though the velocities v,, and Vp (or Vb and Vq) are derived 
from potentials, they appear in convolution with density 
inhomogeneities. For this reason, there is no need for 
explicitly second-order quantities (such as vorticity or 
anisotropic stress) to obtain a nonzero magnetic field, an 
effect similar to the calculation of the kinetic Sunyacv- 
Zel'dovich effect by Ostriker & Vishniac (1986). 

The spectral density. At: k'^ Psik) / {2ttY , obtained by 
numerically integrating 9q and 5q in equations (22) and 
(24), provides both a measure of the magnetic field 
strength on a given scale k and a measure of the energy 
stored in magnetic fields, ps — I-BP/Stt. The results for 
the spectral density of the magnetic field on comoving 
scales from lO^-^Mpc^^ to lO^Mpc"^ at the epoch of 
recombination are shown in Figure 3. The peak of the 
density corresponds to a typical magnetic field strength 
of 10^25_]^Q-24 Q Qn a comoving scale of 0.1 Mpc~^. 
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5. DISCUSSION 

The major result of this paper has been to demon- 
strate that seed magnetic fields of cosmologically inter- 
esting strengths and scales are necessarily generated by 
the same processes that cause structure formation. As 
overdensities in the early universe slowly grow during the 
radiation era, the differing photon interactions with pro- 
tons and electrons create charge separations and currents 
of small magnitudes on all scales. These charge separa- 
tions and currents grow in magnitude as the universe 
evolves, causing magnetic fields to grow as well. The 
magnetic power in a given mode peaks at approximately 
the time of horizon crossing, oscillating and falling slowly 
in amplitude after that. The net result is that, at the 
epoch of recombination (and hence prior to any field 
amplification due to gravitational collapse or dynamo ef- 
fects), seed magnetic fields of magnitude O{10~^'^ G) are 
created by the simple dynamics of charged particles. 

Other recent discussions of the origins of large-scale 
magnetic fields have mentioned the possibility that the 
evolution of cosmological perturbations could be respon- 
sible for their generation (Matarrese et al. 2005; Taka- 
hashi et al. 2005; Gopal & Sethi 2005; Ichiki et al. 
2006). These groups have focused on the effects of 
vorticity, which arises at second order in cosmological 
perturbation theory. In addition to the charge separa- 
tions and currents induced by photon scattering and the 
Coulomb interaction, as discussed in this work, Taka- 
hashi et al. (2005) also considered the contribution of 
photon anisotropic stress to magnetic fields, but found 
that it is uninteresting on large scales. The overall 
strength and scale of the magnetic fields reported by 
these authors seem to differ by large factors from one 
another. The fields derived in Takahashi et al. (2005) 
and Ichiki et al. (2006) appear to be unreasonably large, 
and peak at scales where Silk damping ought to be a 
dominant effect (Matarrese et al. 2005). As shown in the 
subsequent paragraph, their value for a magnetic field 
strength of B r^ 10"^^ G on ©(Mpc) scales is difficuh to 
reconcile with other estimates. It also appears that the 
fields derived in Gopal & Sethi (2005) are insignificantly 
small, as their analysis of vorticity may not include the 
dominant terms in the creation of magnetic fields. For 
reasons discussed in the next paragraph, it appears that 
the estimates for magnetic fields generated by vorticity 
are most accurately given by Matarrese et al. (2005). 
When their results are considered on ~ 10 Mpc scales, 
the fields generated are found to be of physically reason- 
able strengths similar to our own, of order ~ lO"^^ Q. 
Their spectrum is similar to the spectrum of Figure 3 as 
well, although their spectrum peaks at a slightly higher 
amplitude, is narrower in shape, and reaches its maxi- 
mum at a smaller scale. 



The approximate magnitude of the field strengths aris- 
ing from the full calculations of Section 4 can be un- 
derstood intuitively. From the quasi-equilibrium solu- 
tion in equation (22), the fractional charge asymmetry 
(Sq) is of order 10"^'^ at the scale of the peak of the 
power spectrum. When this information is combined 
with the fact that the streaming velocity at decoupling is 
v/c ~ 10^''', a current that yields magnetic field strengths 
of 0(10^^^ G) is produced. The field strengths expected 
from vorticity are somewhat less intuitive. Although for 
an ideal fluid a rotational velocity is not generated from 
an initially longitudinal velocity if there is no vorticity at 
linear order (Hwang 1993), this does not remain the case 
if the imperfect nature of the fluid is properly taken into 
account (Vishniac 1982). A reasonable "upper-bound" 
estimate of the magnetic field that arises from vortic- 
ity can be calculated from the angular momentum of a 
protogalaxy, yielding an angular velocity (to) and there- 
fore a magnetic field. Harrison (1970) showed that mag- 
netic fields of magnitude B ^ 10~^wG would arise in 
the radiation era due to vorticity. From Peebles (1969), 
an estimate of the angular momentum in the baryonic 
component of galaxies is ^ 10^" gcm^ s""'^, which yields 
a magnetic field of _B ^ 10^^^ G for a protogalaxy at 
recombination. This estimate of the field strength is rea- 
sonable as an upper bound, as the angular momentum 
from Peebles (1969) is for the entire protogalactic system, 
but the associated magnetic field energy may appear dis- 
tributed across many smaller-scale modes. Furthermore, 
this estimate is comparable to and consistent with the 
field strengths found in Matarrese et al. (2005). It is 
therefore our assessment that their paper is the most cor- 
rect of the works detailing magnetic fields arising from 
the vorticity generated by cosmological perturbations. 

It should not be surprising that the calculational 
frameworks of this paper and those involving magnetic 
fields arising from vorticity can lead to comparable re- 
sults, as both calculations rely on second order quanti- 
ties. The calculation of a vorticity in the context of cos- 
mological perturbation theory requires a velocity term 
that is explicitly second order, whereas Section 4 of this 
paper requires the convolution of two first order quanti- 
ties. On dimensional grounds, the two methods should 
yield comparable results. The results in this paper rely 
only on fields known to exist, and whether the contribu- 
tion from charge separations and currents or from vortic- 
ity dominates can vary both with scale and from system 
to system. The usual evolution of cosmological pertur- 
bations leads in a straightforward manner to charge sep- 
arations and currents, which necessarily arise from the 
differing interactions between protons and electrons with 
photons. As has been shown in Section 4, magnetic fields 
follow directly from these sources. 

A seed field of 0(10^^'^ G), as generated via the charge 
separations and currents induced by structure forma- 
tion, would provide an excellent candidate for the ori- 
gins of cosmological magnetic fields. While field am- 
plification due to gravitational collapse is negligible at 
the epoch of recombination, this will not be the case at 
all times. At recombination, the universe has only been 
matter-dominated for a brief time, and thus density per- 
turbations have only grown by a small amount in that 
time, leading to an insignificant amplification of the field 
strength. As magnetic flux is frozen in, nonlinear col- 



lapse causes |i?| to increase by many orders of magnitude 
(Lesch & Chiba 1995; Howard & Kulsrud 1997). How- 
ever, the major source of amplification of an initial seed 
field comes from dynamo effects, as discussed in Section 
2. The key to solving the puzzle of the origin of cos- 
mic magnetic fields lies in determining whether the seed 
fields produced by a given mechanism can be success- 
fully amplified into the 0(/xG) fields observed today. A 
major problem with many of the mechanisms that pro- 
duce seed fields is that they produce weak fields at times 
insufficiently early for dynamo amplification to produce 
fields as large as ~ /iG. The Biermann mechanism, for 
instance, can produce seed fields of order ~ 10~^^ G, 
but only at a redshift of z ~ 20. Although those initial 
fields are larger than the ~ lO"^** G fields produced by 
the growth of cosmic structure, the fact that magnetic 
fields from structure formation are in place at z ~ 10^ 
makes them an extremely attractive candidate for the 
seeds of cosmic magnetic fields. As argued by Davis et 
al. (1999), a seed field as small as 10^"^° G at recombi- 
nation could possibly be amplified into a /iG field today. 
Clearly, more work on understanding dynamo amplifica- 
tion is necessary before a definitive solution to the puzzle 
of cosmic magnetic fields can emerge. 

One interesting mechanism worth investigating further 
is for the cosmic seed fields generated by density pertur- 
bations to seed supermassive black holes. It is known 
that the magnetic field energy in active galactic nuclei 
and quasars is comparable to the magnetic field energy 
in an entire galaxy. However, these structures cannot 
generate their own magnetic fields from nothing; they re- 
quire a pre-existing seed field. It therefore appears to be 
a reasonable possibility that the seed fields generated by 
cosmic structure formation could provide the necessary 
fields to seed supermassive black holes. The resultant 
amplification via collapse and dynamo effects could ex- 
plain the origin of large-scale magnetic structures in the 
universe. 

If large-scale magnetic fields exist at the epoch of re- 
combination, they may be detectable by upcoming exper- 
iments. The results shown in Figure 3 provide a predic- 
tion of large-scale magnetic fields at the epoch of the cos- 
mic microwave background. Sufficiently large magnetic 
fields on large scales at recombination may be detectable 
by PLANCK (Lewis 2004; Delabrouille 2004), although 
current estimates of their sensitivity indicate that the 
field strengths predicted in this paper (^ 10^^^ G) would 
be significantly out of range of PLANCK 's capabilities 
(^ 10^^" G). Nonetheless, a knowledge of the field 
strengths at recombination allow for predictions of CMB 
photon polarizations and Faraday rotation, both of which 
may be, at least in principle, observable. 

It is also of interest to note that any primordial charge 
asymmetry or large-scale currents (and therefore mag- 
netic fields) created in the very early universe will be 
driven away by these dynamics. Equation (20) has an 
approximate solution for Sq which is critically (exponen- 
tially) damped at redshift z ~ 10^, capable of reducing 
any pre-existing charge or current by as much as a factor 

of ~ 10"^*^ . As this factor is extraordinarily large, it is 
easy to conclude that any initial 6q or 9q will be driven 
quickly to the value given by equations (22) and (24) at 
the epoch of critical damping. This result is indepen- 



dent of k and ought to be applicable even to a charge 
asymmetry on large scales, perhaps approaching that of 
the present horizon. As a result of this reduction of pre- 
existing charges or currents, the late-time (z < 10^) solu- 
tions obtained in this paper for charge separations, cur- 
rents, and magnetic fields ought to be independent of the 
initial conditions for (5^, 9^, and |_B| in the early universe. 
Overall, the dynamics of ions, electrons, and photons 
during the radiation era necessarily leads to charge sepa- 
rations and currents on all scales, which in turn generate 
magnetic fields. These fields supersede any pre-existing 
fields and are in place prior to substantial gravitational 
collapse. Thus, the creation of charge separations and 
currents from the evolution of cosmological perturbations 



emerges as a promising and well-motivated new candi- 
date to explain the origins of cosmic magnetic fields. 
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